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We report the first observation of tlie baryonic flavor-changing neutral current decay A^ —^ Afi'^ fi~ 
with 24 signal events and a statistical significance of 5.8 Gaussian standard deviations. This mea- 
surement uses a pp collisions data sample corresponding to 6.8 fb~^ at ^/s = 1.96 TeV collected 
by the CDF II detector at the Tevatron collider. The total and differential branching ratios for 
A^ -> Afi+fj.- are measured. We find ^(Ag -> A^+^~) = [1.73 ± 0.42(stat) ± 0.55(syst)] x 10"^ 
We also report the first measurement of the differential branching ratio of B^ — >■ (jjjj,'^ fi~ using 49 



signal events. In addition, we report branching ratios for B"*" 
B ^ K * (892) fi+fj,- decays. 



K+ti+^j.-, B" -> Ky+^i- and 



PACS numbers: 13.25 Hw, 13.20 He, 13.30 -a 
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Rare decays of hadrons containing bottom quarks 
through the process b -^ sii'^fj,", where 6 is a bottom 
quark and s is a strange quark, occur in the standard 
model (SM) with 0(10-^) branching ratios [l|,i|. The b 
and s quarks carry the same charge but different flavor, so 
this process is a flavor-changing neutral-current (FCNC) 
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decay. FCNC decays are suppressed at tree level in the 
SM, and must occur through higher order, and more sup- 
pressed, loop amplitudes. Their suppressed nature and 
clean experimental signature, along with reliable theo- 
retical predictions for their rates [ll [sj, |j], make them 
excellent search channels for new physics. With multi- 
body final states, these decays offer sensitivity to new 
physics in a number of kinematic distributions in addi- 
tion to the total branching ratio. In this Letter we report 
measurements of the total branching ratios of FCNC de- 
cays, as well as their differential branching ratios as a 
function of q^ = A/^ c^, where M^^ is the dimuon in- 
variant mass. Exclusive decays of -B — >■ K^*^ fj,'^ fj,~ have 
been observed by BABAR ji]. Belle 0], and CDF 0. The 
CDF experiment also recently reported the observation 
of 5° — > 0(lO2O)/i"*"/i~ [7|. No significant departure from 
the SM has been found thus far. 

In addition, the study of the baryonic b — > sfi^ fi^ de- 
cays is very important, since the baryonic FCNC decays 
are sensitive to the helicity structure of effective Hamil- 
tonian which is lost in the hadronization of the mesonic 
decays [8j. Although the theoretical calculations of the 
exclusive baryonic b — )■ sjjl^ ii~ decays have large uncer- 
tainties compared to the mesonic decays due to addi- 
tional degrees of freedom in the baryon bound states, the 
measurements of the total and the differential branching 
ratios can help the improvement of the theoretical treat- 
ments. One can also compare the measurements of the 
mesonic b — > s/u+/x~ decays with the baryonic decays, 
which follow the common quark transition. Measure- 
ments of both mesonic and baryonic FCNC decays there- 
fore provide additional tests of the SM and its extensions. 
However, no b baryon FCNC decay has been observed 
and there are few experimental constraints on their decay 
rates. The A° — > Au +/x~ decay is considered promising 
in this respect |8l4ll| and experimentally accessible since 
the branching ratio is predicted as (4.0 ± 1.2) x 10~^ [lOj. 

The data sample used in the measurements reported 
in this Letter corresponds to an integrated luminosity of 
6.8 fb~^ from pp collisions at a center-of-mass energy of 
y/s = 1.96 TeV collected with the CDF II detector be- 
tween March 2002 and June 2010. The K^ -^ A^l+^l-' 
decay is reconstructed and measurements are made of 
the total branching ratio and the differential branching 
ratio as a function of g^. Besides the updated branch- 
ing ratios of B^ — > (t)fi~^fi~, 5+ -^ K~^^~^iJ.~, and 
B° — >• K*{892)'^^^fi~, we report the branching ratios 
of B° -^ K°^i+fi- and B+ -^ K*{892)+n+j2- which are 
measured for the first time in hadron collisions. We also 
report the first measurement of the differential branch- 
ing ratio as a function of q^ of B^ — > 0/i+/i^. To cancel 
the dominant systematic uncertainties, decay rates for 
each rare channel Hb — >■ h^'^ ii~ are measured relative 
to the corresponding resonant channel Hb — > J /iph with 
J /if) — > fX^ fjT , used as a normalization, where Hb repre- 
sents the b hadron and h stands for A, 0, K^ , Kg, K*^, 



and K*'^. Charge-conjugation is implied throughout the 
Letter. 

The reconstruction of the exclusive b hadron events 
starts with a dimuon sample selected by the online trig- 
ger system [l^ of the CDF II detector 13 1. The trigger 
system utilizes information from muon detectors and the 
central outer tracker [ij]. Muon chambers CMU and 
CMX [ll cover |r?| < 0.6 and 0.6 < \tj\ < 1.0, respec- 
tively, [16J. The CMP muon chamber covers |?7| < 0.6 
and is located behind the CMU and an additional steel 
absorber. The dimuon trigger requires a pair of oppo- 
sitely charged particles with a momentum transverse to 
the beam line pr > 1.5GeV/c, which are matched to 
track segments in the CMU or CMX chambers. At least 
one of the muon tracks is required to have a CMU track 
segment. The trigger also requires that the dimuon pair 
satisfies either L^y > 100 fj,m, where the transverse decay 
length Lxy is the flight distance between the dimuon ver- 
tex and the event primary vertex JlT], or px > 3.0 GeV/c 
and matched segments in both CMU and CMP chambers 
for one of the muon candidates. 

Offline event selection starts with the triggered dimuon 
pairs. Each offline track is required to satisfy more strin- 
gent requirements on the number of hits used to recon- 
struct the track. The dimuon selection requirements used 
in the trigger are repeated with the higher quality of- 
fline tracks. The decay length and invariant mass of 
each dimuon pair are calculated after a vertex fit using 
the muon tracks. Dimuon pairs are classified according 
to their invariant mass M^^. Dimuons from FCNC b 
hadron decays are required to be inconsistent with de- 
caying from J/ip {i/j{2S)) mesons by requiring q^ values 
outside the window of 8.68 (12.86) < g^ < io.09 (14.18) 
GeV^/c^ [7\. The J/tp candidates are required to have 
Mfj,fj, within 50 MeV/c^ of the known J/ip mass [l^l • 

The A° —i' A/i+/i^ candidates are selected by combin- 
ing the dimuon pairs with A baryons reconstructed from 
decays A — >■ pir^ . The p7T~ pairs are required to have 
invariant mass consistent with the known A mass jl8| . 
Pt ^ 1-0 GeV/c, and a vertex displaced from the dimuon 
vertex. The transverse momentum of the A^ candidate is 
required to be greater than 4.0 GeV/c. Candidates with 
an invariant mass calculated from two or three daugh- 
ter particles compatible with J/tp, tp{2S), D°, £)+, D+, 
or Ac masses are rejected to remove backgrounds from 
these charm- hadron decays [7|. The B^ — > ipiX^ jjT can- 
didates are reconstructed from dimuons together with a 
pair of oppositely-charge kaons consistent with a (j) de- 
cay with a selection similar to that of A° — >■ A/j,"'"/i~. 
The i3°'+ — >• /C°'+/i+/i^ candidates, where KP^^ is one of 
{K^ ,K^,K*'^ ,K*^}, are formed from a dimuon com- 
bined with up to three charged tracks. The K'g me- 
son is reconstructed in its 7r"'"7r~ final state by requir- 
ing the dipion mass to be consistent with the known 



Kl 



mass 



18[. Details about the reconstruction of the 
decays of K*° -^ K+i:' and (j) -^ K'^K' can be 



Cross-feed between A° 



A/z"'"/x and 



found in Ref. 0- 

B° — )■ Kgfj,'^iJ,~ is suppressed by evaluating the momen- 
tum imbalance of A and Kg daughters [19| . We utilize 
the correlation between invariant mass and the asymme- 
try a = {ql - ql) / {ql + ql), where g^ ^"^ is the longitu- 
dinal momentum of the positive (negative) decay product 
relative to the direction of the A or K^. We reject can- 
didates satisfy -0.26 < -l.QM{Kl) + \a\ < -0.15 for 
Ksf^+fi~ and 4.73 < 3.6M(A) + \a\ < 4.78 for Afi+ ^i' . 
We remove 76 (90)% of the cross-feed while the signal loss 
is 11 (7)% for A^+/i^ {Kgfi'^^~). A residual cross-feed 
contamination of 0.1% (0.6%) to the A/x+^~ {Kgfi~^^~) 
signal is considered as a systematic uncertainty. To fur- 
ther optimize the event selection, an artificial neural net- 
work (NN) classifier is trained using simulated signal 
events and background events taken from Hb mass side- 
bands (0.1-0.36 GeV/c^ far from the known Hj, mass) in 
data. Some kinematical distributions of the simulated 
signal, e.g., the transverse momentum of b hadron, and 
the energy depositions of muon candidates in the elec- 
tromagnetic and hadron calorimeters, are corrected us- 
ing scale factors extracted by comparing simulation to 
data in the normalization channels. We use 70% of the 
sideband events for the training, and use the remaining 
events to check that the NN does not bias or over suppress 
the mass distribution. The optimized NN threshold is de- 
termined to maximize the average expected significance 
of the branching ratio, using many kinematic observables 
including transverse momentum, invariant mass, vertex 
fit qualities and muon identification qualities [Tj. 

The signal yield of the A° — > A/i+/i^ candidates is ob- 
tained by an unbinned maximum likelihood fit to the 
AjJ invariant mass distribution with the signal proba- 
bility density function (PDF) parametrized by Gaus- 
sian distributions using simulated signals and the back- 
ground PDF modeled by a linear function. We fix the 
A^ mass width for the rare decay while it is fioated 
for the normalization channel. Different mass width be- 
tween data and the simulated signal is corrected by mea- 
sured mass width ratio of the normalization channel be- 
tween data and the simulated signal. The signal region 
is defined within ±40 MeV/c^ from the world average Aj] 
mass [7\. The statistical significance is obtained through 
a likelihood-ratio test between the signal plus background 
and background-only hypotheses interpreted assuming it 
distributed as a x^ variable. The invariant mass dis- 
tribution of the Aj] — )■ A^"*"^" candidates is shown in 
Fig- HI In the signal region, we observe 24 ±5 events from 
A" — >■ A/j,+^~ decays while the total number of the signal 
candidates is 34. The statistical significance of the signal 
s corresponds to 5.8 Gaussian standard deviations. The 
signal yields of B'^ — > (pfX^ fi" and other FCNC B me- 
son decays are obtained by a similar procedure as that 
of Aj — > Kfi^ fi^ . Each channel uses independent NN 
weight and PDF. The fit range for B+ and B^ decays 
is from 5.18 to 5.70GeV/c^ to avoid the region of 5.0- 
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FIG. 1: Invariant mass of (a) A° -5> A/i+^", (b) B° -5> 
(ji^+M", (c) B+ -^ K+fi+fi-, (d) B° -^ K*°^j.+fi-, (e) B° -^ 
Ksfi+fj.~ , (f) B+ -^ K'+fi+fi" , with fit results overlaid. The 
histograms are the data. Solid, dashed-dotted, and dotted 
curves show the total fit, the signal PDF and the background 
PDF, respectively. 



5.18GeV/c2, 
ground from 
contribution 
to the muon 
between B'^ 
ulation, and 
of S^ ^ (j)^i- 
shown in Fig 



which is dominated by the feed-down back- 
multibody decays of b hadrons. While the 
from charmless Hh decays is negligible due 

identification, we estimate a 1% crosstalk 
-> K*^^^fi~ and B^ — >• 0^+//~ using sim- 
correct for it. Invariant mass distributions 
^fi^ and other FCNC B meson decays are 

[1] and signal yields are listed in Table H] 



TABLE I: Summary of observed yields, the statistical signif- 
icance s, and the relative efficiency Erci. 



Mode 


Nh^+u- 


s{a) 


Nj/^h 


£rcl 


a;; ^ A/i+M" 


24 ±5 


5.8 


1740 ± 50 


0.33 ±0.01 


5° ^ 0/.+/.- 


49 ±7 


9.0 


4560 ± 80 


0.56 ±0.01 


B+ -^ K+fi+ij.- 


234 ± 19 


13.7 


72200 ± 300 


0.41 ±0.01 


BO ^ K'°,j.+fi- 


164 ±15 


13.7 


28300 ± 200 


0.45 ± 0.02 


B° ^ K°s^^+f^- 


28 ±9 


3.5 


9470 ± 90 


0.47 ±0.01 


B+ -^ K*+ii+ii- 


20 ±6 


3.5 


4560 ± 80 


0.38 ±0.02 



The branching ratios of A° — >• A/z^/x", B^ — >• (/)/i+/z~, 
and B — )■ K^*'fi~^iJ,~ are calculated by comparing their 
signal event yield to that of the normalization decay 
modes Ag -^ J/7/-A, S^ ^ J/i/>(/), and B -^ J/'ipK'^*\ 
where J/tp — >■ M^/^ i after the reconstruction efficiency 
correction: 



<r- 2 






N, 



hfi^ fj, 



B{J/^^|l+^l-) 



N 



J/iph 



frcl 



(1) 



where Ni^^+^- is the hfi^fi'' yield, Nj/^j^ is the 
J/^h yield for the normalization channel, and Erci = 
^hfi+fj." /^j/tlih is the relative reconstruction efhciency de- 
termined from the simulation. The calculated relative 
and absolute branching ratios are listed in Table [TTl The 
absolute branching ratios are obtained using world aver- 
ages of the J/tph decay rates 



of B" -^ K^fi+^i- and B+ 



m 



K*+n+fi- 



The branching ratios 
are measured 
for the first time in hadron collisions. 

The dominant sources of systematic uncertainty are 
the scale-factor reweighting of the simulated signal (the 
trigger efficiency near the threshold) which ranges from 
0.5% to 4.0% (0.8% to 7.2%), depending on the channel. 
We estimate the former uncertainty from the comparison 
of the relative efficiencies with and without reweighting 
and the latter uncertainty from the different px require- 



ments for each trigger. 



In the A° 



-^ A^+/i case we 



consider an additional uncertainty of 6.6% due to the 
unknown A'^ — > J/ipA polarization. 

For the absolute branching ratio measurements we as- 
sign the uncertainties on the world average B(Hh — >■ 
J/iph) [18| or the most recent measurement [20[ . Contri- 
butions from other sources (e.g., background PDF shape 
or the decay model of the simulated event) are minor 
(0.3%-3.4%). 

The combined branching ratio is calculated by assum- 
ing isospin symmetry and using the B^ and B^ total 
widths [l8|. These numbers are consistent with our pre- 
vious results [7| , B-factory measurements [5|, |6[ , and the- 
oretical expectations [3, [lO| . 

We also measure differential branching ratios with re- 



TABLE II: Measured branching ratios of rare modes. First 
(second) uncertainty is statistical (systematic). The last two 
values are for the isospin average. 



Mode 


Relative Z3(10~^) 


Absolute Z3(10~'') 


a;; ^ A^+M" 


2.45 ±0.59 ±0.29 


1.73 ±0.42 ±0.55 


B« -^ <^M+M" 


1.13 ±0.19 ±0.07 


1.47 ±0.24 ±0.46 


B+ ^ K+^i+^i- 


0.46 ± 0.04 ± 0.02 


0.46 ± 0.04 ± 0.02 


B« ^ i^*V+M" 


0.77 ± 0.08 ± 0.03 


1.02 ±0.10 ±0.06 


bo ^ K°^i+^i' 


0.37 ±0.12 ±0.02 


0.32 ±0.10 ±0.02 


B+ -^ K*+^l+fl~ 


0.67 ±0.22 ±0.04 


0.95 ±0.32 ±0.08 


B -^ Kfj,+^~ 


- 


0.42 ± 0.04 ± 0.02 


B ^ K'fi+fi- 


- 


1.01 ±0.10 ±0.05 
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FIG. 2: Differential branching ratios of (a) A5 — >• A/i+/i , (b) 
B" -> 0M^M", (c) B+ -^ K+fi+fi-, and (d) B° -^ A'*V^At"- 
The points are the fit result. The solid curves are the SM 
expectation [ll. I lOl. |21| . The dashed line in the A^ — >■ Afi^ fj,~ 
plot is the SM prediction normalized to our total branching 
ratio measurement. The hatched regions are the charmonium 
veto regions. 



spect to q^ . We divide the signal region into six bins in q^ . 
We fit the signal yield in each q^ bin. In each fit, we fix 
the mean of the Ht mass and the background slope to the 
value from the global fit, so that only the signal fraction 
is allowed to vary in the fit. Figure [2] shows the differ- 
ential branching ratios for A" -^ A/i+/i^, S^ -^ (f)fi'^ij.~, 
and B —>■ K^*' ^^ fi^ . For illustration, we superimpose 
the SM expectations, which are based on the formula in 
Ref. [ij), with the form factors in Ref. 21 1, 
the case of A° — >■ A/i^/i" 
The cusp at q"^ 



except for 
decays which follows Ref. jlO| . 
7GeV^/c^ is due to a change in pa- 



rameter approximations. Tables IIIII and IIVI summarize 
the differential branching ratio measurements. The two 
bottom rows in each table show the results for the semi- 
inclusive bins which are included with ranges covering 
theoretically well-controlled regions. 

In summary, we have updated our previous analysis of 
the flavor-changing neutral current decays b — > s^,^ pT 
using data corresponding to an integrated luminosity of 
6.8 fb~ and adding new decay channels. We report the 
first observation of Ajj — > A/i+/i~ and measure the total 
and differential branching ratios of this decay with re- 
spect to q^ . We also measure the total and differential 



branching ratios of B — > K^*^ 



fi^fi 



and B° -^ 0^-1 



Ai 



with respect to q-. All measurements are consistent and 
competitive with other results, and the differential mea- 



TABLE III: DifTerential branching ratios of A" — > A/i'^ /i , _B+ — > _ftr+/i+/i , B" — > K'^fi^/i , combined B — > J^/i^/i , in units 
of 10^ . The Qmax IS 20.30 (23.00) GeV /c for Afi'^ fi^ (/C/i^/i^). The first (second) uncertainty is statistical (systematic). 



g" (Gcvyc^) as; ^ A^+^ 



B^ 



K 



j^nn 



B 



TT 



Ky+fi- 



B -» Kfi^fi' 



[0.00, 2.00) 

[2.00, 4.30) 

[4.30, 8.68) 

[10.09, 12.86) 

[14.18, 16.00) 

[16.00, 



2 
I ymax/ 



0.15 ±2.01 ±0.05 
1.84 ±1.66 ±0.59 
-0.20 ±1.64 ±0.08 
2.97 ±1.47 ±0.95 
0.96 ±0.73 ±0.31 
6.97 ±1.88 ±2.23 



0.36 ±0.11 ±0.03 
0.80 ±0.15 ±0.05 
1.18 ±0.19 ±0.09 
0.68 ±0.12 ±0.05 
0.53 ±0.10 ±0.03 
0.48 ±0.11 ±0.03 



0.312 ±0.372 ±0.024 0.33 ± 0.10 ± 0.02 

0.929 ± 0.485 ± 0.070 0.77 ± 0.14 ± 0.05 

0.663 ±0.510 ±0.052 1.05 ± 0.17 ± 0.07 

-0.030 ± 0.223 ± 0.005 0.48 ± 0.10 ± 0.03 

0.726 ±0.257 ±0.055 0.52 ± 0.09 ± 0.03 

0.214 ±0.182 ±0.016 0.38 ± 0.09 ± 0.02 



[0.00,4.30) 2.65 ±2.52 ±0.85 1.13 ± 0.19 ± 0.08 1.268 ± 0.622 ± 0.096 1.07 ± 0.17 ± 0.07 
[1.00,6.00) 1.27 ±2.08 ±0.41 1.41 ± 0.20 ± 0.10 0.980 ± 0.614 ± 0.076 1.29 ± 0.18 ± 0.08 



TABLE IV: Differential branching ratios of B° -s- (pfi'^ fi' , B° -> K'^fi+fi', B^ 



in units of 10 
(systematic). 



The g^a^ is 18.90 (19.30) GeV'^/c^ for <^m^M" {K* h'^ 



— > K*'^ fi^ jj. , and combined B — >■ K*jj.'^fi , 
The first (second) uncertainty is statistical 



q' (GeVVc^) 


B^' -^ <^m"^^" 


B" ^ 7^*V+M- 


B+ -^ K'+fi+ii- 


B -^ K*fi+n~^ 


[0.00, 2.00) 


2.78 ±0.95 ±0.89 


1.80 ±0.36 ±0.11 


1.30 ±0.98 ±0.14 


1.73 ±0.33 ±0.10 


[2.00,4.30) 


0.58 ±0.55 ±0.19 


0.84 ±0.28 ±0.06 


0.71 ±1.00 ±0.15 


0.82 ±0.26 ±0.06 


[4.30, 8.68) 


1.34 ±0.83 ±0.43 


1.73 ±0.43 ±0.15 


1.71 ±1.58 ±0.49 


1.72 ±0.41 ±0.14 


[10.09, 12.86) 


2.98 ±0.95 ±0.95 


1.77 ±0.36 ±0.12 


1.97 ±0.99 ±0.22 


1.77 ±0.34 ±0.11 


[14.18, 16.00) 


1.86 ±0.66 ±0.59 


1.34 ±0.26 ±0.08 


0.52 ±0.61 ±0.09 


1.21 ±0.24 ±0.07 


[16.00, gLx) 


2.32 ± 0.76 ± 0.74 


0.97 ±0.26 ±0.07 


1.57 ±0.96 ±0.17 


0.88 ±0.22 ±0.05 


[0.00,4.30) 


3.30 ±1.09 ±1.05 


2.60 ±0.45 ±0.17 


2.01 ±1.39 ±0.27 


2.53 ±0.43 ±0.15 


[1.00,6.00) 


1.14 ±0.79 ±0.36 


1.42 ±0.41 ±0.12 


2.57 ±1.61 ±0.40 


1.48 ±0.39 ±0.12 



surements of i?^ — > (pfi^n^ and A|J -^ A/i+/i^ arc the 
first such measurements. At present there is no evidence 
of discrepancy from the SM prediction. 
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